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tonian H -I- 9{t)AH and into correlation functions. With the help of recent work on correlation 
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Witten diagram in the two-dimensional AdS-Vaidya background, and find agreement with the 
conformal quantum mechanics result. 
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1 Introduction and summary 


Partially motivated by the AdS/CFT correspondence, there has been recent progress in under¬ 
standing correlation functions in conformally invariant quantum mechanics (CQM), relevant for 
the AdS 2 /CFTi case. In this case, the 50(2, l)-isometry of the bulk AdS 2 manifold manifests as 
a 50(2,1) conformal invariance of a quantum mechanical theory on the boundary. 

The benchmark model of conformal quantum mechanics, viewed as the 0 = 1 limit of 50(0 -|- 
1,1) conformally invariant scalar field theories with a potential in O spacetime dimen¬ 

sions, was studied by de Alfaro, Fubini and Fur Ian (dAFF) [T]. The action is invariant under 
50(2,1) i.e. 5L(2,M) transformations, with the generators satisfying, after quantization, the com¬ 
mutation relations 

i[D,H]=H, i[D,K] = -K, i[K,H]=2D (1.1) 


of a so(2,1) Lie algebra. There are many other quantum mechanical realizations of (1.1) besides 
the original dAFF model. For the AdS 2 /CFTi correspondence, we expect the relevant model to 
have many {N) interacting degrees of freedom, with the bulk calculations corresponding to a large 
N, strong coupling limit. However, the full picture of the AdS 2 /CFTi correspondence is not yet 
well understood, due to complications arising from the AdS 2 fragmentation and backreaction when 
the geometry is reached from string theory in a form AdS 2 x X where A is a compact space [2]. 
Recently, these issues were studied in context of 1-1-1 dimensional models of dilaton gravity coupled 
to matter [3]. There one can also construct solutions for an ingoing null matter pulse into a vacuum, 
creating a black hole, so that the two-dimensional spacetime is described by the AdS-Vaidya metric. 

The null collapse to a black hole, or the AdS-Vaidya spacetime has been much studied as an 
analytic holographic model of one type of a global quench in the QFT on the boundary (see [1] for 
an early idea, and e.g. la 0 HI El E] for early papers on AdS-Vaidya and holographic quench). In 
particular, the model enables one to study nonequilibrium Green functions for a strongly coupled 
QFT, as these can be computed from the bulk geometry. However, it is not very transparent how 
exactly the quench is realized on the boundary. Furthermore, especially in higher dimensions the 
Green functions do not usually have a simple analytic form. In this paper, we are interested in 
studying the holographic quench in 0-1-1 dimensions, and the structure of the Green functions. We 
work in a simplified setting, assuming that the bulk spacetime is described by the two-dimensional 
AdS-Vaidya metric. For the dual theory at the boundary, we rely on the recent work cni, where 
the authors studied how to construct correlation functions in conformal quantum mechanics. Using 
only the so(2,1) algebra, without specifying the underlying theory, they constructed operators 
which satisfied some, but not all properties of conformal primaries, and a vacuum state which 
was not fully invariant under (1.1). However, despite these shortcomings when the vacuum and 
the operators are combined to correlation functions, the result satisfies expected transformation 
properties for a conformally invariant theory. 

In addition to the general motivation of extending the study of quenches modeled by AdS-Vaidya 
spacetimes to the AdS 2 /GFTi case, we have specific motivations and goals: 1) it is interesting to 
study to what extent non-equilibrium correlation functions can be computed working directly on 
the boundary theory; 2) we would like to understand in more detail what kind of a quench the 
AdS-Vaidya bulk spacetime corresponds to in the boundary; 3) 3-pt and 4-pt correlation functions 
have been studied at thermal equilibrium m and more recently in the context of semiclassical 
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limits by bulk geodesic methods na Ea [H], but not yet in the context of global quenches. In 
this work we interpret the AdS 2 -Vaidya spacetime to be realized in CQM as a sudden change of 
the Hamiltonian. We compute analytic results for thermal and non-equilibrium 2-pt, 3-pt and 4-pt 
functions in quenched conformal quantum mechanics, by using the results of uni, and the 2-pt and 
3-pt functions also from a bulk AdS 2 -Vaidya background calculation, finding agreement with the 
results from the CQM computation. 

This paper is organized as follows: Section 2 reviews some key concepts of conformal quan¬ 
tum mechanics and various choices of time evolution; in Section 3 we introduce the AdS 2 -Vaidya 
spacetime as a model of a quench, which we find to be realized in CQM as a sudden change of the 
Hamiltonian, and then compute non-equilibrium two-point and 3-point correlation functions both 
by the holographic method from the bulk spacetime and by a direct calculation in CQM. Section 
4 extends the CQM calculation to the 4-point function, and ends with some brief comments. 


2 Review 


We first review some key features of so(2,1) representation theory. The energy eigenvalues (of 
H) are continuous and the eigenstates are all non-normalizable. In addition, there is no state 
that would vanish under all three generators i.e. that would be invariant under all the symmetry 
transformations. Let us then consider another complete set of orthonormal states. Moving to a 
new basis of operators with the commutators 


R = 

2 \ a 


— aH ] 2ziD , 


L, = 

2 V a 

where a is a constant with the dimension of timaT the commutation relations become 


( 2 . 1 ) 


[i?, L±] — ±L±, [L_, L_|_] — 2R , 


( 2 . 2 ) 


so that they form the Cartan-Weyl basis of the so(2,1) algebra. The generator R is compact with 
a discrete set of eigenstates. The discrete lowest weight representation can then be constructed 
using L± as raising and lowering operators. Explicitly: 


R\n) 

rn 

{n'\n) 

L±\n) 

|n) 


rn\n), 

vq + n , n G N 

\Jrn{rn ± 1) - ro{ro - l)|n ± 1), and 


r(2ro) 


re!r(2ro -|- n 


<L+)Q0) 


^The constant a played the role of an infrared regulator in [T]. 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 
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where rg is a parameter, the lowest weight of the representation, also connected to the eigenvalue 
^o(^o ~ 1) of the Casimir invariant C = ^(HK + KH) — = E? — L_|_L_ of the algebraj^ In 

particular, the lowest eigenstate |0) = |n = 0) plays the role of a vacuum, called the i?-vacuum. 
For a generic rg, the eigenvalue equation implies that the ii-vacuum cannot be annihilated by all 
the so(2,1) generators. 

One can also construct a continuous basis of states, |r). Using the conventions in mm. the 
operators take the form 


H = -i 


dr 


D = -i(r—+ro) 

dr 

K = -i(r^^+2ror) 
dr 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


Constructing the operator R in this basis, one can relate the |n) and |r) bases with the help of a 
differential equation: 


{T\R\n) = rn{T\n) = 2 [(a + ay + 2ro^J (r|n) 

Mr) = (r|n) = {-!)” [E 23 ±!!l] ^ (gg)’'" 


’’0 


( 2 . 11 ) 


The constants in the expression are determined by the condition that the raising and lowering 
operators operate in the same way in both bases. 

The |r) basis is not orthonormal, but the overlap of two states has the form 


00 / 

(ti|t 2) = ^n(Ti)/3;(r2) = r(2ro) ( T 

n=0 


2(ri - T 2 ) 


2ro 


( 2 . 12 ) 


the same as for primary operators of conformal weight rg in a conformal theory. In m. the authors 
gave an explicit interpretation of (2.121 as a vacuum 2-point function, constructing an operator 0(t) 
which reproduces the |r) states operating on the R-vacuum, |r) = 0(r)|0). Explicitly, 


0(r) = Ai(r)exp(-a;(r)L+) , 


(2.13) 


with 


so that 


V{r) = [r{2r„)]>/'^(54l^y’'” , 

, , a + iT 

u{t) = . . 

a — IT 

(2.14) 

( 7 - 21 x 1 ) = (0|0^’(r2)0(ri)|0) . 

(2.15) 


^The value of ro depends on the theory. For example, in the inverted harmonic oscillator model introduced in [T], 
ro is related to the value of the dimensionless coupling constant A in the potential term —\q~^. 
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There are some issues with this form. While it does produce the correct state |r), the action 
of on the i?-vacuum just returns the R-vacuum multiplied with the normalization factor N{t). 
This would cause the time-translation invariance to be lost when considering (0|0(Ti)0^^(r2)|0). 
One could consider an alternative definition for O which would still produce the same |r) state^ 
We could use 


0(r) = N{t) exp[—t(;(r)L+] exp[—a;*(r)L_] 


With this definition, 0(r) would produce the same |r) states but this time 0^(r) produces the 
state |r)*, i.e. a state with complex conjugated coefficients. This operator is almost Hermitian, 
apart from the complexity of N(t). This alternative definition is not necessary for the calculations 
in this paper. 

However, there is another alternative form, suggested in m- They suggested 

r) = 

' 2^ro 

as it correctly reproduces the |r) state when acting on |n = 0). The leftmost exponential gives the 
time evolution from the |r = 0) state which the rightmost exponential prepares from the R-vacuum. 


(2.16) 


2.1 Alternative time evolutions 

So far we have chosen H = as the time evolution generator. Alternative choices have been 

considered previously in [laiisiiiz] and this enables one to compute Green functions periodic in 
imaginary time, associated with a finite temperature background. We motivate the alternative time 
evolution generators by reviewing the isometry of two-dimensional anti-de Sitter spacetime AdS 2 . 

Three metrics are often used in the context of AdS 2 : the Poincare metric (the metric written in 
the Poincare coordinate patch covering a part of the spacetime manifold), the global metric (using 
the global coordinates patch that covers the full AdS covering space manifold), and the black hole 
metric (using a coordinate patch, which after periodic identification in imaginary time defines an 
anti-de Sitter black hole) (see [TB] and Figure 3 therein for an illustration of the patches). The 
metrics have the form 


Global : ds'^ 


cos^ 6 


i-dT^ + de^), 


—oo < T < oo, 


TT 

2 


< 9 < 


TT 

2 


(2.17) 


Poincare : ds'^ = -^{—dr'^ + dz‘^), —oo < r < oo, 0 < z < oo 


Black hole : ds^ = 




/ 




du^ 


\ 


1 - 


, —oo < t < oo, 0 < n < uh 


(2.18) 

(2.19) 


^There are some other remaining issues too. It can be shown that exp[R/Ar]|r) = |t -t At), but 
exp[iJlAr]0(r) exp[—R/Ar] 7^ 0 (t -|- At) which makes the Heisenberg picture a bit problematic. In addition, 
the R-vacuum expectation value {0|O(r)|0) does not vanish, unlike what happens in CFT for a primary operator - 
the R-vacuum is not conformally invariant. 
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and the coordinates are related by the transformations 


r ± z 

z 

v{u, t) 


2a tan 




{ T ± e ±-)/2 

2u 


1 + cosh{v/ uh) + {u/uh) sinh(r;/Miif) ’ 


t — Uh artanh ( — ) 

\uhJ 


T — z = 2uh tanh 


2uh) ’ 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


Here, we used different symbols for the variables of different choices of metrics for the sake of clarity. 

The Killing vector fields that generate the SO{2, 1) isometry of AdS 2 are represented in Poincare 
coordinates by 


H = -idr 

D = -i {rdr + zdz) 


K = -i 




dr + 2Tzdz 


(2.23) 


After transforming to global coordinates, the Killing vector that generates time translations in 
global time T is the linear combination (see also m) 


K 

Hq = aH + — = —iOt 
4a 


(2.24) 


Likewise, after transforming to black hole coordinates, time translations in the time t are generated 
by the combination 

Hbh = H- KI{2uh? = -idt . (2.25) 

On the (conformal) boundary of the spacetime, the relation between the different time coordinates 
reduces to 

T = 2a tan(T/2) = 2uH^suih.{t/2uH) (2.26) 

The linear combinations (2.24) and ( |2.25 1 are the alternative time-evolution generators considered 
in m- Moving to conformal quantum mechanics, the generators are represented as infinitesimal 
translation operators when acting on the time basis states, e.g. H\t) = —For the generators 
Hg,Hbh, one needs to define new time states |r), \t) [T5]. Since |r) = O(r)|0), and (in correlation 
functions) 0(t) transforms under coordinate transformation like a primary of weight ro, one defines 
new time states \t) by 


it> = ( 


dtV<^ 
di 


\t = t{i)) , 


or in our case explicitly: 


IT) = a'’o(cos(T/2))-2^o|r = 2atan(T/2)) 

\t) = (cosh(t/(2tt//)))“^'’°|r = 2uh / {2uh))) ■ 


(2.27) 


(2.28) 

(2.29) 
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The different choices of a Hamiltonian are then represented in the different time state basis as 


Hplr) = 

(2.30) 

HclT) = 

(2.31) 

to 

II 

1 

(2.32) 


inherited from the time evolutions in the three customary AdS 2 coordinate patches. The two- 
point function G 2 (r 2 ,ri) = {t 2 \ti) can then be related to two-point functions (T 2 IT 1 ) and {t 2 \ti) in 
terms of a simple rule. Although the operator 0(r) is not a primary one, it transforms as one in 
correlation functions evaluated with respect to the i?-vacuum state |0) = |n = 0), 


(felfi) ^ {0|6t(*2)O{(,)|0) = (0| . 


(2.33) 


This leads to 


(r2|ri) = (o|oJ,(r2)OG(Ti)|o) = 

{t2\tl) = (0|0]j^(t2)0s/f(tl)|0) = 


r(2ro)a^^° 

[4z sin( ^^~^^ )]^’"o 
r(2ro)a2^o 


[4tt//fsinh(^|;^)]2^o 


(2.34) 


3 AdS-Vaidya geometry and holographic quench 

We now move to consider a holographic model of a quench, and incorporate it into conformal 
quantum mechanics. In gauge-gravity duality, a quench in the strongly coupled theory on the 
boundary has a holographic dual interpretation in the bulk AdS geometry. Perhaps the simplest 
and most popular model that has been studied is the AdS-Vaidya spacetime. It describes lightlike 
collapse of matter into a black hole in AdS£)+i space. The holographic dual interpretation is that 
the theory on the D spacetime dimensional boundary is initially in a vacuum state, then is at 
t = 0 instantaneously sourced in a homogeneous manner into an excited state which then time 
evolves into thermal equilibrium. The initial non-equilibrium state is somewhat special, because 
expectation values of all local operators thermalize instantaneously. The non-equilibrium nature 
is only revealed by considering expectation values of non-local operators, such as 2-point functions 
or Wilson loops. The two-dimensional AdS-Vaidya spacetime is even more special, because the 
boundary has no space directions. In this paper, we are intereted in correlation functions. 2- 
point (autocorrelation) functions in AdS 2 -Vaidya background have been computed in [2D1 El] . We 
rederive the result by a simpler calculation using a coordinate transformation. We will then show 
that it is very simple to obtain in conformal quantum mechanics. 

The AdS 2 -Vaidya spacetime (the Penrose diagram is depicted in Figure see also |3], Figure 
1, for a related spacetime illustration) is specified by the metric 


ds^ = 


^2 L 


2 

— 9{v)^^dv‘^ — 2dvdz 


(3.1) 
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Figure 1: Penrose diagram of the AdS 2 -Vaidya spacetime. 


For u < 0, the metric (3.1) is the metric of AdS 2 written in terms of a lightcone time coordinate 


V = T — z. For u > 0, the metric (3.1) is the metric of an AdS 2 black hole spacetime, with 
dv = dt — dz/{l — {z/zh)"^) (we used u instead of z in (2.19)). It is useful to note that the black 


hole spacetime metric can also be transformed to the AdS 2 metric with a change of coordinates 
(that also acts on the boundary) 


2z 


z = 


1 + cosh.{v/ zh) + {zj zh) s\nh.{v/ zh) ’ 


V = 2zH^9uVih{v/{2zH))- 


(3.2) 


We can also combine the vacuum and black hole regions of the Vaidya spacetime by a continuous 
piecewise coordinate transformation 


z = 


2 , 

2z 


l+cosh{v/ zh)+{z/Zh) sinh(ti/2/i-) ’ 


u < 0 
u > 0, 


and 


V = 


V, u < 0 

2zH^Wih.{v/{2 zh)), u > 0 , 


so that the AdS 2 -Vaidya metric reduces globally to that of AdS 2 , 

ds^ = ^ ~ “^dvdz^. 


(3.3) 


(3.4) 


(3.5) 


For correlation function calculations in the AdS-Vaidya background, it is simple to work with the 
{z,v) coordinates. Final results are then obtained by performing the inverse coordinate transfor¬ 
mation in the end. 








3.1 AdS-Vaidya 2-point function 

The (near boundary) vacuum two point function in AdS 2 is well known. Let us consider for 
simplicity a massless bulk scalar field, 


G^p\zi,Vi] Z2,V2) ~ - 


ZlZ2 


-IT —{V2 — Vl)‘^ + 


(3.6) 


We then perform the inverse coordinate transformation to the original coordinates v, z and obtain 
the in-in (vacuum) two point function in the AdS 2 -Vaidya. Note that the bulk Klein-Gordon 
equation of motion uniquely fixes the full two point function with the initial condition that it must 


for hi < 0 and h 2 < 0 agree with the vacuum two point function. Clearly this is the case for (3.6). 

Let us then consider the region h 2 > 0 and hi < 0. Then, we use the inverse coordinate 
transformation, and take the boundary limit with the scaling prefactors [ 1 ] to find the boundary 
two point function for an operator with scale dimension A = 1. 


t2 


G''p\{t 2 ,ti) = lim (^ 12 : 2 ) ^^^(-Zi, ti; 22 /cosh^ 22 ;//tanh 


2:1,22 —^0 

1 1 


2zh ' 


t2 , 


1 


cosh^ ^ vr —{ti — I tanh + ie 
1 1 
TT —(ti cosh ^ — I sinh ' 


(3.7) 


where we use the notation I = 2zh- From now on, we will be using the parameter I which gives us 
the temperature T = l/( 7 r/). This result was obtained earlier in j 2 ()ll 21 ll^ using more complicated 
methods. For V 2 > 0 and vi > 0, we simply obtain the thermal two point function 


G%{t2,ti) = 


1 


1 


(3.8) 


vr —(I sinh L ^)2 _|_ 

Similarly, for an operator with a generic scale dimension A 7 ^ 1, the exponent in the denominator 


changes, e.g. (3.6) becomes 


G%{t2,ti) = - -^-- 

^ [“(^1 cosh ^ — I sinh y)^ + ie] 


(3.9) 


3.2 


Holographic quench and correlation functions in conformal quantum me¬ 
chanics 


The AdS-Vaidya spacetime corresponds to a quench in the boundary theory, but its precise re¬ 


alization has not been transparent. From the discussion in Section (2.1) we learn that the two- 


dimensional AdS-Vaidya bulk geometry corresponds to a sudden change of the time evolution of 
the system, thus in conformal quantum mechanics the quench corresponds to the Hamiltonian 


H = Hp + e{t)AH , 


(3.10) 
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where Hp generates time evolution with respect to the Poincare time r and A.H = —K/{2uh)‘^, 
so that after the quench the Hamiltonian becomes Hbh, generating time evolution with respect to 
the black hole time coordinate t. The two time coordinates have a common origin r = t = 0 at the 
quench. 

Further, in 1+1 dimensional conformal field theory with a quench, it is non-trivial to compute 
non-equilibrium correlation functions [23]. On the other hand, in simple free quantum mechanical 
systems such as the harmonic oscillator, one usually works in the Heisenberg picture where in- 
in vacuum correlation functions can be calculated by working out the appropriate Bogoliubov 
transformation. In conformal quantum mechanics, at least for the holographic quench of interest 
here, the situation is simpler since we do not need the explicit form of the Hamiltonian of the 
system which would be needed to evaluate the Bogoliubov transformation. We only need to apply 
the rule how the operators (or more specifically, correlation functions) transform under a conformal 
transformation. For example, consider the non-equilibrium two-point function with t 2 > 0, the 
operator 0{t2) is inserted after the quench, and ti < 0, insertion before the quench. For the 
insertions, we defin^ 


n /'+ - J , for t > 0 

UVaidya^t) j , for t = T < 0 


(3.12) 


and recall that r = tanh(t//) for t > 0, r = t for t < 0. We then use the operator Ovaidyait) for 
the two insertions. Thus, in order to calculate the correlation functions, we only need to apply 
the rule (2.33) how the operators (in fact, correlation functions) transform under a conformal 
transformation, to arrive at 


{^\Ovaidyai^‘^)Ovaidya{tl)\0) - F(2ro) 


a 


2ro 


2i[ti cosii (^) — I sinh (^)] 


(3.13) 


The result is in agreement (after adjusting the overall normalization) with the above bulk calculation 
(in the limit e = 0), when we match the lowest weight with the scale dimension, ro = A. 


3.3 Three-point functions 


In conformal quantum mechanics, it is equally straightforward to compute the three- and four- 
point functions in the holographic quench background. This is a new result, since previous studies 
of holographic quenches have been focusing on two-point functions. We consider the three-point 
functions which have already been calculated in the zero temperature conformal quantum mechanics 

nnnE]: 


{t2\(I){t)\ti) = (0|0(r2)(()(T)0(ri)|0) 


A 

(n - T2)^^0 “'^(t - ri)'^(r - T2)^ ’ 


(3.14) 


“^Our justification for this definition is the following. Recall that the operator 0{t) is defined through the state 
\t) = Cl(t)|0). We are lead to solve the non-equilibrium quench problem by solving the time dependent Schrodinger 
equation 

-idt\t) = H\t), (3.11) 


where H is the quenched Hamiltonian (3.10l. In the Appendix, we show that (3.111 is solved by \t) = Ovaidya (t)|0). 
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where <p is a primary operator with scale dimension 6 and ^ is a constant 


A = (0|<^(0)|0) 



2ro+S 

r(2ro)a2"°. 


(3.15) 


(Note that (0|i;i)(0)|0) needs not vanish, as the R-vacuum is not invariant under all conformal 
transformations.) 

Making a coordinate transformation, t = Zartanh (j), and using the conformal transformation 
rules in the correlator, the three-point function becomes 


{t2\(l){t)\ti) sinh^'-o-^ (^) ’ 


(3.16) 


which is the expected finite temperature result. 

Now, we do a quench at t = 0, e.g. we turn on the temperature at t = 0, and compute the 
three-point function with t 2 > 0, < 0. First, we consider t < 0. In this case. 


{t2\(l){t)\ti) 


A 


{t — tiY (t — l tanh (y))”^ (ti — I tanh (y))^^° ^ cosh^'’“(^) 

A 

{t — ti)^ (tcosh(^) — /sinh (y))^ (ticosh(^) — Zsinh (y))^^° ^ 


(3.17) 

(3.18) 


In the denominator of the last expression, we see the familiar factors from the thermalizing two- 
point function. 

If t, t 2 > 0, < 0, the 3-point function becomes 


{t2\(j){t)\ti) = (3.19) 

_ _ A _ 

(/tanh (j) — ti)*^ (/tanh [j) — / tanh (y))”^ (ti — / tanh (y))^^° ‘^cosh^^°(^) cosh^‘^(j) 

^ _ f _ 

(/sinh(|) — /i cosh(^))‘^sinh'^ (cosh(^)ti — / sinh(^))2'’o-<5 

We will next compare these results against a holographic derivation of 3-point functions from the 
AdS 2 -Vaidya background. We begin by reviewing some facts of the calculation in an AdS 2 vacuum 
background ~ since we are interested in in-in vacuum correlation functions, we need to adopt a 
Keldysh contour method, which we present next. (Note: we focus only on the leading contribution 
to the 3-point function, ignoring the additional contributions associated with backreaction which 
could be computed if bulk 1-|-1 (dilaton) gravitational dynamics would be included as in [3].) 


3.4 AdS 2 vacuum 3-point functions with Keldysh contour 

As in our previous bulk analysis, to keep matters simple we consider a massless scalar field. For 
the three-point function we add self-interaction terms and start from the action 

5 = + ...). (3.20) 


11 










X2 

X 


Xi 

—X— 


Xs 

—X— 


-<- 

Figure 2: The complex time (Keldysh) contour for the in-in 3-point function. 


The bulk to bulk Feynman propagator of cj) is given by 


G^P{xi,X 2 ) = - —log 


-{t2 - Ti)^ {zi - Z2f + ie 
-{t 2 - Ti)^ + {zi + Z 2 y + ie 
^1^2 , ^^, 3 


- , 2 - + 0{z^), (3.21) 

TT-{t 2 - ny + zj + le 

where we denote xj = {zj, tj). In the end we are interested in the limit 2:2 0 and thus, we will 

work with the order 0 ( 2 ^ 2 ) term above. Also, we will in the following need the Wightman two- 
point function, which can be obtained from the Feynman one using the identity {4>{xi)(p{x2)) = 
0{t2 - Ti){T{(t>{xi)4>{x2))) + 0{ti - T2){T{4){xi)4>{x2)))* , and is given by 


G^l\xi,X2) = - 


Z 1 Z 2 


TT -{t2 - ri)2 + z\ + i{T2 - Tl)e 


+ 0 (^ 2 )- 


(3.22) 


The 3-point function can be now calculated using perturbation theory in A. The calculation would 
be simplest in Euclidean time, but since we will later consider a non-equilibrium situation, which 
is inherently real time, we will show how to calculate the vacuum 3-point function in the bulk real 
time formalism. The boundary correlator is obtained as a limit of the bulk correlator by using the 
extrapolate dictionary. 

The bulk 3-point function, defined in terms of the Heisenberg picture field operators 0//, can 
be written in terms of the Dirac/Interaction picture operators cj), using a complex time contour as 

G’P{xi,X2,X3) = {T{(I)h{xi)(I)h{x2)4>h{x3))) = (Tc ((^(xi) (()(X2) (()(x3)e“*^ -/c ^ 

(3.23) 

where Tq denotes time ordering along the complex time (Keldysh) contour shown in Fig. 


Expanding (3.231 to first order in A, and using Wick’s theorem we obtain 


G^^(xi,X 2,X3) =-iA [ (fx{Tc{(l){x)(j){xi))){Tc{(l){x)(l){x2))){Tc{(l){x)cj){x3))). (3.24) 

Jc 

Writing out the contour integral gives 

G^^(xi,X 2,X3) =-iA J (fx{T{(f){x)(j){xi))){T{(l){x)(j){x2))){T{(l){x)(j){x3))) (3.25) 

+ i^J d^x{(j){x)(j){xi)){(f){x)(j){x2)){(f)ix)(j){x3)), (3.26) 

where now the time integrals run from r = — 00 to r = -|-oo. 

Thus, the 3-point function is a sum of two Feynman diagrams of the form shown in Fig. 
with the lines denoting in the first case time ordered two-point functions, and in the second case. 
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Figure 3: The first order Feynman diagram contributing to the scalar three-point function. 

Wightman two-point functions. Using the known two-point functions gives 

3 


gP{xi,X 2,X3) = -i\ J (fx^J\g\f[ 


1 




TT -(r - Ti 


•U z'^ + ie 


-h 


iX [ (fx^\ 


ZiZ 


*=i ^ + Z^ + i(T - Ti)e 

The boundary three-point function in the "extrapolate" dictionary is given by 

G^a(U,t2,t3) = lira {ziZ2Z3)~^G^^'>{xi, X 2 , xs). 

Thus, we finally obtain the bulk 3-point function as the sum of the following integrals 

3 


(3.27) 


(3.28) 


G^p^{n,T2,T3) = -iX I dzdTzYl 


1 


1 


+ iX dzdt z 


TT -(r - ri)2 z"^ X- ie 

3 1 


1 


(3.29) 


vr -(r - ri)2 + z"^ + i{T - rj)e 
We will first compute the r integrals, from the residues at the poles of the integrands. The Wight¬ 


man two-point functions in the second line of (3.29) have poles at 


T = Ti ± z + ie- 


(3.30) 


As the poles are all located in the upper part of the complex r plane, we can close the integral 
contour from below without encountering any poles. As the integrand vanishes as at large r, 
the integral vanishes. Thus, the 3-point function reduces to the contribution from the time ordered 


2-point functions, in the first line of (3.29). This is the expected result in the vacuum state, where 


the in-in and in-out formalisms are expected to agree. The time ordered two-point functions have 
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poles at 


r = r, ± z ± ie. 


(3.31) 


Performing the r integral by closing the integral contour from the upper half complex r plane gives 




/O 


(3.32) 


where Tjk = Tj — Tk- The z integrals are now elementary and can be performed using the identity 

/■“ dz 1 


/o {z + a){z + b) a-b 


'<)■ 


(3.33) 


leading to 


G%{ti,T2,T3) = 


A 1 


log ( 


27 r 2 T12T13T23 

Using log(—1 — i 6 ) = —i7rsign((5) for (5^0, gives finally 

A3) ' . . A 1 


_ 1 _ .^ Ji + r| + r| - T 1 T 2 - T 1 T 3 - T 2 T 3 


'ri2n37'23 


G'V,a(u,T2,r3) = -i 


2vr |ri2r23Ti3| 


(3.34) 


(3.35) 


3.5 Vaidya 3-point functions 


Next, consider the 3-point function in the AdS-Vaidya spacetime. Using the real time Schwinger- 
Keldysh formalism, with the above two-point functions gives the same integral expression as the 


ground state (3.27), but now in the barred coordinates. One difference to the vacuum case is that 


the integral over light-cone time now has an upper limit at v = 2 zh due to the coordinate relation 


(3.4). It is easy to see that the integration region can be continued all the way to u = 00 without 


changing the value of the integral, as in this region, the integrands in the first and second lines of 


(3.27) identically cancel each other. This follows from unitarity, as the Schwinger-Keldysh contour 


can be extended forwards in time without changing the result for the correlation function. Thus, 
the same calculation of the integral in the vacuum 3-point function goes through and we obtain 
the near boundary 3-point function 


^(3) , A 2:12:223 

'^F ~ I - - 

2 tt \T12T23T13 


(3.36) 


Now depending on whether fj is before or after the collapse, we get different results. When all of 
the points are in the region fj <0, we obtain the vacuum result. For fi < 0, f 2 < 0 and fs > 0, we 
obtain using the extrapolate dictionary 



■ ^ ^ 1 
27r cosh^ X |(ti — t 2 ){ti — I tanh x)(t 2 — I tanh ^)| 


(3.37) 
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On the other hand if two of the points ^2 and are located after the shell, we obtain 


27r cosh^ ^ I 

And finally when all of the points are located after 

\ 1 

- X 


tne 


^(3) . 
'-'F,d - 


X 


A 1 

27r cosh^ ^ cosh^ ^ cosh^ ^ 

1 


|(/ tanh Y 


-1 tanh 


^)(/tanh ^ 


/tanh Y)(/tanh y — / tanh y)! 


(3.38) 


(3.39) 


4 Four-point functions in finite temperature and with a quench 


In conformal quantum mechanics, it is almost as straightforward to derive the non-equilibrium four- 
point functions as the three-point functions. For completeness, and for possible future reference, 
we end with this calculation. (For example, we expect that the structure of similar non-equilibrium 
3- and 4-point functions in higher dimensional field theories reflects those of conformal quantum 
mechanics and can be reduced to them in an appropriate (equal space) limit. However, we expect 
the computation to be more involved in higher dimensions.) At zero temperature, for fields (j) and 
(p with dimensions 6 and 5, respectively, the four-point function can be evaluated as (note, the 
corresponding expression in m contains a typo) 




F(2ro)x’-°2Fi(^,J;2ro;x)(0|(()(0)|0)(0|<^(0)|0) 
2S+S+2ro (t24)'^-^0 (il2)^+'’° (/34)'^+'’“ ’ 


(4.1) 


where Uj = U — tj and 


^12^34 

X =-, 

^13^24 

is the conformally invariant ratio familiar from the usual conformal field theories. In the above, we 
have set the scaling parameter a to unity. Here, the x dependent factors remind us of the model 
dependent functions of 4-point functions in conformal field theories. It is noteworthy that the 
above expression can be obtained by using only a single conformal block. Now, we do the change 
of variables to the thermal coordinates. 





sinh( ) sinh( ^ 


(4.3) 


and 


(0|<^(0)|0)(0|(^(0)|0)F(2ro)x[° 2i^i(^,^; 
sinh(^)^“^o sinh(^)'^“’'o sinh(^)'^+^o sinh(^)'5+^o sinh(^)2^0“^“'^ 

Now for the real thing, thermalization. Now, we set the quench moment at / = 0. Also, t 4 > 0 
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and ti < 0. We consider three different cases. First, t 2 < 0 ,t 3 < 0: 




( 0 |</.( 0 )| 0 )( 0 |<^( 0 )| 0 ) 

(tl3)^-’'0(ti2)'5+^0 


(4.5) 


_r(2ro)x^° 2Fii5,5; 2ro; _ 

(t 2 sinh(^) — Z sinh(^))‘^“^o(t 3 cosh(^) — / sinh(^))^+^o(ti cosh(^) — / sinh(^))2^0“‘^“'^ 


with 


Then, t 2 < 0, ^3 > 0, 


^ ti 2 (cosh(^)t 3 - /sinh(^)) 

ti3(cosh(^)t2 - ^sinh(^)) 


(4.6) 


{tl\4>{t2)^{t3)\t4) 


( 0 |,/.( 0 )| 0 )( 0 |<^( 0 )| 0 ) 

(ti2)<5+^o sinh(^)^+^o 


(4.7) 


_r(2ro)x^o 2Fi{6, ~6-, 2ro; _ 

(ticosh(^) — Z sinh(^))'^“^o(t 2 sinh(^) — Z sinh(^))^“'’o(ti cosh(^) — Z sinh(^))2^o-<5-5 


with 

_ ZZi2sinh(^) 

^ (cosh(^)ti — Zsinh(^))(cosh(^)t 2 — Zsinh(^)) 
Finally, when both t 2 ,t^ > 0, we have 


(4.8) 


where 


(tl|(/>(t2)(^(t3)|t4) 


_ (O|0(O)|O)(O|.^(O)|O) __ 

(ti cosh(^) — Z sinh(^))^+’'o(Zi cosh(^) — Z sinh(^))^“'’o 


F(2ro)x^° 2 Ti(^, 2ro; x,)2-^-^-^^°Z-^-^ 
(ticosh(^) — Z sinh(^))2^0“‘^“'^ sinh(^)^+^o sinh(^)‘^“^o ’ 


(4.9) 


_ (ti cosh(^) — Zsinh(^)) sinh(^) 
^ (ti cosh(^) — Zsinh(^)) sinh(^) 


(4.10) 


The results are of course what one expects, the structure of the four (and three)-point functions 
change in a systematic way at different points of time with respect to the quench moment. 

For further study, it would be interesting to include the backreaction corrections to the 3-point 
functions [3], or to study how the 3- and 4-point functions in CQM are recovered from higher 
dimensions. The non-equilibrium 2-point function agrees with the equal space limit of the 2-point 
function in 1-|-1 dimensions, in the limit of large operator dimension when the 2 -point function is 
well approximated by the geodesic approach. The analytical results for the correlation functions 
would be interesting to compare with the equal space limits of cZ-|-1 correlation functions for cZ > 1, 
because the blackening factor in the bulk black hole metric is dimension dependent. It would also be 
interesting to see how the (semiclassical) limit of 4-point functions computed in higher dimensions 
[lainiiii] reduce to the 4-point function in CQM. 
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Appendix 


We motivate the definition Ovaidya{t) in (3.121 by checking that the time evolution of the state 
|r) = Op(r) jo) connects to that of \t) = OBp(t)| 0 ) after the quench. 

Recall that 


= Op(r)| 0 ) = e 


_ irHp aHp 


|o) 


(4.11) 


where the second exponential prepares the state jr = 0 ) and the first factor continues the 

time evolution (also to earlier times). The state jr = 0) = |t = 0) at the quench. After the quench, 
the time evolution continues with the new time translation generator Hpp, so it should evaluate 
to the state 


\t)' = Opp(t)| 0 ) = , 


(4.12) 


On the other hand, we had defined the state \t) to be 



ro 


T = r(t)) 


[cosh(t//)] jr = / tanh(t//)) 


(4.13) 


where I = 2uh- Thus, we need to show that the two states are the same, {t)' = \t). To establish 
that, we study how they time evolve. On one hand. 


.2ro 


-idt\t) = 

= cosh“^(t//) 


sinh(t//)[cosh(t//)] ^ + cosh 


\t) 


27*0 

Hp + cosh(t//) sinh(t/Z) 


It). 


(4.14) 


On the other hand. 


- idt\t)' = HbhW = {Hp + ^K)\t)' , 


(4.15) 


where 7 = We then show that the above two first order equations are in fact identical. 

Since the initial state |t = r = 0) = e“^r|0) is the same, that ensures that |t)' = |t). 

Starting from 


K\t) 


-i [r'^dr + 2 rod 1 ^) = {j^Hp - i 2 ror^ |t) 
— fcosh“^(t//) — 1 ) Hp + -— p tanh(t/Z) 

_7 \ / 7 < 


It) 


(4.16) 
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we obtain 


iHp + jK)\t) = 


cosh ^(t//)i7p + tanh(t//) 


\t) = -idt\t) 


so that \t) indeed satisfies the same first order equation as \t)'. 


(4.17) 
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